Resonant inelastic x-ray scattering study of holon-antiholon continuum in SrCuC>2 
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We report a resonant inelastic x-ray scattering study of charge excitations in the quasi-one- 
dimensional Mott insulator SrCuC^. We observe a continuum of low-energy excitations, the onset of 
which exhibits a small dispersion of ~ 0.4 eV. Within this continuum, a highly dispersive feature with 
a large sinusoidal dispersion (~ 1.1 eV) is observed. We have also measured the optical conductivity, 
and studied the dynamic response of the extended Hubbard model with realistic parameters, using 
a dynamical density-matrix renormalization group method. In contrast to earlier work, we do not 
find a long-lived exciton, but rather these results suggest that the excitation spectrum comprises a 
holon-antiholon continuum together with a broad resonance. 

PACS numbers: 78.70.Ck, 78.30.-j, 71.10.Fd, 75.10.Pq 



The separation of spin and charge degrees of freedom 
is one of the most important and fascinating properties of 
electrons in strongly correlated systems in one dimension. 
In particular, it is well known that in the one-dimensional 
(ID) Hubbard model, the low-energy physics is domi- 
nated by collective excitations of decoupled charge and 
spin degrees of freedom called holons and spinons, re- 
spectively 0. Experimentally, if one creates a hole by 
removing an electron, this hole is expected to decay into 
a spinon and a holon, which can be studied with an- 
gle resolved photoemission spectroscopy (ARPES) 0- 
The situation is different for so-called "particle-hole" 
probes, such as optical spectroscopy, resonant inelastic 
x-ray scattering (RIXS), and electron energy loss spec- 
troscopy (EELS). In these experiments, the total charge 
is conserved in the scattering process, so that an electron 
is simply moved from one site to another, creating a hole 
and a doubly occupied site. The decay of the hole cre- 
ates a holon and a spinon, while the double-occupancy 
decays into an antiholon and a spinon. Since photons 
and electrons strongly couple to the charge sector, the 
behavior of holon-antiholon pairs can be studied with 
these particle-hole probes 
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The so-called corner-sharing chain cuprates S^CuOs 
and SrCu02 are both charge-transfer insulators; that 
is, they have insulating gaps of ~ 2 eV arising from 
strong electron correlations. Since their crystal structure 
is highly anisotropic, the electronic structure remains ID 
over a wide temperature range. Only at a very low tem- 
perature does magnetic order set in, due to the small in- 
terchain coupling (T N w 2 K for SrCu0 2 @). Based on 
the commonly-used measure of quasi-one-dimensionality, 
Tjv/J ~ 10~ 3 , these compounds can be regarded as 
among the best realizations of quasi- ID systems [f|. The 



quasi-lD nature of both SrCu02 and S^CuOs has been 
studied extensively with various experimental techniques, 
including magnetic susceptibility pf, ARPES 2|, Raman 
scattering ||, optical spectroscopy 0, EELS ^j, RIXS 
0, and neutron scattering [j"fT |. 

In this Letter, we report a detailed study of the mo- 
mentum dependence of the low-energy charge excitations 
in SrCu02, utilizing the RIXS technique. We observe 
a continuum of excitations arising from the creation of 
particle and hole pairs. Within this continuum, a well- 
defined spectral feature with a large sinusoidal disper- 
sion (~ 1.1 eV) is observed. We have also measured the 
optical conductivity cr(u>) and carried out a dynamical 
density-matrix renormalization group (DDMRG) calcu- 
lation of the extended Hubbard model. A comparison of 
the RIXS spectra with both the theoretical results and 
cr(u>) yields a consistent picture of charge excitations in 
this material: We find two dispersive spectral features, 
in contrast to the earlier RIXS study 4J. The first fea- 
ture is the onset energy of the holon-antiholon contin- 
uum, which exhibits small dispersion (~ 0.4 eV). The 
second spectral feature is a more dispersive broad reso- 
nance that is the remnant of the strong-coupling exciton. 
This result contrasts earlier EELS work, in which the ex- 
citation spectrum was interpreted as being dominated by 
a bound exciton at large momenta 

In RIXS experiments, the incident x-ray energy is 
tuned to near the absorption edge of the particular ele- 
ment of interest, so that the inelastic scattering intensity 
of certain electronic excitations is resonantly enhanced. 
In addition, this gives an element specificity that is valu- 
able in studying complex materials such as the cuprate 
compounds H El El El E3 The RIXS experiments 
were carried out at the Advanced Photon Source on the 
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FIG. 1: Pseudocolor plot of RIXS intensity as a function of 
momentum and energy transfer (Ei — 8982 eV). 



undulator beamline 9IDB with the same setup used in 
previous studies ^ij. The scattering plane was vertical 
and the polarization of the incident x-ray was kept per- 
pendicular to the Cu-0 plaquette. Note that the corner- 
sharing CuC>2 chain runs along the c-direction. A single 
crystal of SrCuC>2 was grown using the traveling solvent 
floating zone method. The crystal was cleaved along the 
(0 1 0) plane and mounted on an aluminum sample holder 
at room temperature in an evacuated chamber. 

In order to find the resonance condition, we first car- 
ried out a study of the incident energy (Ei) dependence of 
the RIXS intensity. There are two features that show res- 
onance behavior in the energy-loss (lu = Ei — Ef) spectra, 
at oj ~ 3 eV and u> ~ 6 eV, respectively. The intensity of 
the 3 eV feature shows large enhancements at Ei « 8982 
eV and Ei s» 8990 eV, while the 6 eV feature is reso- 
nantly enhanced at E t w 8990 eV and E 4 « 8996 eV. 
We note that these three resonance energies correspond 
to the peaks in x-ray absorption spectra near the Cu K- 
edge |l5j. The intensity at 6 eV likely arises from both 
the Cu K/?5 emission and a charge-transfer type excita- 
tion commonly found in the RIXS studies of cuprates 
[ill fl2]| . Detailed studies of the independence will be 
reported elsewhere, and here we fix Ei = 8982 eV, and 
focus on the momentum dependence of the 3 eV feature. 

Our main results are shown in Fig. ^ in which the 
RIXS intensity is plotted as a function of momentum 
and energy transfers. The pseudocolor scale of the ob- 
served intensity is shown on the right side in units of 
counts per second. The momentum transfer is shown in 
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FIG. 2: The dispersion relation of the peak position (filled 
symbol) and the onset energy (open symbol) of the RIXS 
spectra. The solid line is ui(q) = 3.07 — 0.55 cos(g). Calcu- 
lated peaks and onset energies in dynamical density- density 
correlation function, N(q,ui), are also plotted as plus and X 
symbols, respectively. The dashed line denotes the spinon 
dispersion relation as described in the text. 



reduced units of q/2n along the chain direction, and cor- 
responds to the (0 11 q) position. Thus <7/27r = —0.5 
is the ID Brillouin zone boundary, while q/2n = is 
the zone center. The lower cut-off of u> = 0.5 eV shown 
in Fig. ^ is due to the large quasi-elastic contribution, 
and u) = 6 eV is the upper limit of accessible energy 
transfer in this setup. The plotted intensity is the raw 
RIXS intensity without any absorption corrections. Al- 
most constant factors were obtained from three different 
normalization methods |lfj| . enabling us to compare the 
absolute intensities between different momentum trans- 
fers directly, except for q/2ir = — 1. The anomalously 
low overall intensity at q/2n = —1 arises from the very 
low angle of the incoming photons at this q. 

The most prominent feature of the RIXS spectra shown 
in Fig. is the highly dispersive feature at ui m 3 eV. 
The momentum dependence of the peak position of this 
feature, shown in Fig. [21 exhibits a clear sinusoidal dis- 
persion with a bandwidth of 1.1 eV. However, closer in- 
spection of Fig^reveals the presence of additional spec- 
tral weight on the low-energy side around q/2n = —0.5. 
This low-energy intensity around the zone boundary is 
evident for the q = ir data in Fig. |3 where the scat- 
tered intensity is plotted on a logarithmic scale. Since 
the momentum resolution at this position is about 0.1 
reciprocal lattice units, it is unlikely that the additional 
low-energy intensity is due to a resolution effect. Thus, 
the best description of the observed spectra is that of a 
continuum of excitations, in which the dispersive feature 
resides. Note that this interpretation, particularly near 
q = 7T, is in contrast to previous EELS and RIXS results. 
In their EELS study of S^CuOs, Neudert et al. reported 
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that the sharp spectral feature observed near q = ir is a 
bound exciton mode, due to a strong inter-site Coulomb 
interaction The earlier RIXS study by Hasan et al. 
attributed the dispersive feature to the Mott gap edge; 
that is, a lower bound of particle- hole continuum 4]. Our 
RIXS results suggest that the continuum starts at much 
lower energy, and that the particle-hole pairs do not form 
a bound exciton at q = n. The onset energy of the spec- 
tral weight extracted from each scan is shown in Fig. |2| 
as an arrow. The dispersion of these onset energies are 
plotted in Fig. [21 Note that the bandwidth of the onset 
energy dispersion (~ 0.4 eV) is much smaller than that 
of the peak position, as discussed later. 

In concert with the RIXS data, we have also studied 
the optical properties of SrCu02, measuring the reflectiv- 
ity of a cleaved surface between 6 meV and 6 eV, with the 
polarization along the c-axis. Optical conductivity, cr(u>) 
calculated from the reflectivity data using the Kramers- 
Kronig relation with an anchoring process, is plotted as 
a solid line in Fig. 0] and compared with the q = RIXS 
data. The validity of the Kramers-Kronig calculation 
was checked by directly measuring the complex dielectric 
function with an ellipsometry method. In the frequency 
range of interest (1 ~ 4 eV), identical spectra were ob- 
tained from both techniques. The q = RIXS data is ex- 
pected to be different from that of a(tu) — aTm[e(0, w)], 
since the response function measured in RIXS is believed 
to be close to a dielectric loss function, Im[— l/e(q, ui)}. 
Instead, one should focus on the onset energy of the spec- 
tral weight, which should be independent of both the 
particular response function being measured and matrix- 
element effects. Both the RIXS data and the cr(ui) main 
feature exhibit onset energies around 1.7 eV, which cor- 
responds to the continuum edge at the zone center. 

To understand the observed continuum and the dis- 
persive spectral feature, we have studied a simple one- 
band half-filled ID extended Hubbard model (EHM) 
[l7L Hil 0, , in which the nearest-neighbor repulsion 
V is explicitly included, in addition to the usual t and U. 
It is clear that a full theoretical description of SrCu02 re- 
quires a multi-band model based on the Cu 3d and O 2p 
orbitals. However, a precise and reliable determination of 
dynamical correlation functions in such models is still a 
formidable challenge. Several previous studies suggested 
that much of the essential physics of both Sr 2 Cu03 and 
SrCuO? at low energies can be captured by the EHM 
Ullll E3 We have therefore used a DDMRG 
algorithm 20] to determine the optical conductivity and 
density response function in this model for realistic val- 
ues of U and V on chains of up to 128 sites. All previous 
studies either dealt with unrealistic parameter values 
or suffered from severe finite-size effects caused by the 
very small system sizes amenable to exact diagonaliza- 
tion methods [2l|. 

The parameters t, U and V of the EHM were fixed by 
comparing DDMRG calculations of the optical conduc- 
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FIG. 3: Comparison of the RIXS spectra (circles) and Af(q, w) 
(solid lines) at various momenta. To facilitate a direct com- 
parison, Af(q, u>) is calculated with an experimental broaden- 
ing of r) = 0.5t and divided by sin 2 (g/2). Arrows denote the 
onset of the RIXS spectral weight. 



tivity <t(ui) to the experimental data and requiring com- 
patibility with the value J « 0.23 eV for the Heisenberg 
exchange coupling J obtained from neutron scattering 
[To). As shown in Fig. 0] an excellent fit for the low- 
energy peak in (t(ui) is obtained for t — 0.435 eV, V/t — 
1.3, U/t — 7.8. For these values we obtain J « 0.24 eV 
from both 1/U expansions and DDMRG computations of 
the spinon dispersion in the single-particle spectral func- 
tion. Note that the calculated <r(ui) in Fig. ^agrees with 
the measured spectra on an absolute scale. 

The excitation spectrum of the EHM with these pa- 
rameters consists of scattering states of gapped, spinless 
charge =Fe collective modes - (anti)holons, and gapless, 
chargeless spin ±1/2 collective modes - spinons. In con- 
trast to the strong-coupling results reported in Ref. ^3 > 
no excitonic excitations are found in any part of the 
Brillouin zone. The current operator is charge neutral 
and hence only excitations with equal numbers of holons 
and antiholons contribute to cr(u>). Most of the spectral 
weight of the low-energy peak in <j(lu) is attributed to ex- 
cited states with one holon and one antiholon each. The 
contributions of the spin degrees of freedom to cr(u)) is 
small [l9j. 

The RIXS data is more difficult to model as at present 
there is no completely satisfactory theory of RIXS. How- 
ever, it is clear that the RIXS process involves charge 
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FIG. 4: The optical conductivity a(ui) at T = 90 K is plotted 
as a solid line, and compared with RIXS spectra obtained at 
the zone center. The dashed line is the DDMRG calculation 
of cr{uS) with an experimental broadening of n = O.lt. 



fluctuations that may be represented by the low energy 
excitations of the EHM dEHHI. We have therefore de- 
termined the dynamical density-density correlation func- 
tion Af(q, uj) by DDMRG. The results are shown in Fig.|3] 
as solid lines. At small momenta Af(q, uj) is proportional 
to q 2 and thus we use a normalization factor of sin 2 (g/2) 
to compare A/*(g, uj) at various momenta. We note that on 
a lattice the Fourier transform V(q) of the Coulomb po- 
tential is approximately given by l/sin 2 (g/2). Therefore, 
the renormalized Af(q, u>) in Fig. [31 is probably close to 
the dielectric loss function Im[— l/e(q, uj)} oc V(q)Af(q,uj) 
of the EHM. At q = there is a broad "band" of states 
with a peak around 2.3 eV. As we move towards the 
zone boundary, this feature narrows significantly. Im- 
portantly, however, even at q = u the peak has an in- 
trinsic width, that is, it is not a long-lived exciton mode. 
Furthermore, there is sizable intensity below the peak at 
q = 7r. This is in contrast to the strong-coupling U — > oo 
limit |17| . where jV(q ~ is dominated by an exci- 

tonic holon-antiholon bound state. Our interpretation is 
that as U is reduced to smaller, more realistic values, the 
exciton acquires a finite lifetime and turns into a broad 
holon-antiholon resonance. 

Let us try to relate these findings to the RIXS data. In 
Fig- El the normalized Af(q, to) is directly compared with 
the RIXS spectra at equivalent g-positions. The onset 
of the RIXS spectrum shown in Fig. seems to lie in 
the area where Af(q, uj) begins to show appreciable spec- 
tral weight. The onset energy of jV(q,ui) is plotted in 
Fig. [3 as a function of q (crosses). The onset dispersion 
of j\f(q,ui) follows the dashed line, which is the single- 
spinon dispersion relation shifted by a constant energy: 
§ J | sin(g)| + 1.9. This spinon-like dispersion of the on- 
set energy of A/*(g, uj) is consistent with the low-energy 
field theory prediction [22|. Qualitative features of this 



spinon-like dispersion of the onset energy of Af(q, uj) ap- 
pears to be consistent with features of the RIXS data 
(open circles, Fig.|2J), such as the bandwidth of ~ ttJ/2, 
and the q = ir point being a local minimum. In addition, 
we have compared the dispersion of the peak positions. 
As is indicated in Fig. |3 the peak dispersion in Af(q, uj) 
is very similar to that of the RIXS data. However, such a 
quantitative agreement between the RIXS peak positions 
and Af(q, uj) may well be coincidental, and further calcu- 
lation of the full RIXS response function is necessary to 
draw any firm conclusions. 

In summary, we have studied the energy and mo- 
mentum dependence of low-energy charge excitations in 
quasi-lD SrCu02, using both optical spectroscopy and 
resonant inelastic x-ray scattering. We observe a con- 
tinuum of excitations which we associate with the cre- 
ation of holon and antiholon pairs. The onset of the 
holon-antiholon continuum exhibits spinon-like disper- 
sion. Within this continuum, a well-defined spectral fea- 
ture with a large sinusoidal dispersion is also observed. 
By comparing the observed spectra with a DDMRG cal- 
culation of ID EHM, we interpret this well-defined fea- 
ture as a broad holon-antiholon resonance due to the 
inter-site Coulomb interaction. In contrast to earlier 
studies, we find that the holon-antiholon pairs do not 
form an exciton, and the continuum starts at much lower 
energy than previously reported. 
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